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ABSTRACT 


Measurements  of  the  differential  light  scattering  cross- 
section  in  the  ocean  show  that,  while  the  volumetric  differential 
scattering  cross  section  d0/d(2  is  highly  peaked  in  the  forward 
direction,  the  mean-square  scattering  angle 

<  ©2  >  -  (eT)“1  /  o2(dg/dn)d£i 

is  not  small.  Consequently,  the  many  previous  theoretical 

o 

treatments  based  on  the  smallness  of  <  0  >  are  not  applicable. 
The  differential  scattering  cross-section  is  well  approximated  by 
the  form 


0(0)  -  0Teo[ir20(02  +  02)]"1 


where  0^  *  .12  radians  (7°)»  A  solution  to  the  radiative  transfer 
equation  is  found  for  this  differential  cross  section.  Many  useful 
results  can  be  obtained  analytically;  e.g.,  the  spatial  spreading  of 
an  initial  pencil  beam  is  given  by 


T 

Ar  -  0Q  /  dt'[z(x)  -  z(t')]  , 


where  Ar  is  the  full  width  at  half  maximum,  and  z(t)  Is  the 

z 

depth  associated  with  optical  depth  t  ■  /  0T(z')dz'  .  This  rate 

of  spreading  exceeds  that  of  treatments  b2sed  on  a  small  mean-square 
scattering  angle. 
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LIST  OP  ILLUSTRATIONS 


Figure  1 


Figure  2 


Measurements  of  volumetric  differential  scattering 
cross-section  and  the  theoretical  fit,  formula  (2) 
[From  Fig.  15,  N.  G.  Jerlov,  Marine  Optics, 
Elsevier,  (1976),  p.  38}... . 

Geometry  relevant  to  degradation  of  a  point-source 
image.  The  origin  of  z  is  at  the  point  source. 
The  vector  %  is  the  component  of  a  unit  vector 

A 

n  which  is  perpendicular  to  the  z-axls.  The 
imaging  system  is  far  removed  and  records  only 
radiation  with  &  -  0  . 


INTRODUCTION 


The  propagation  of  light  in  clear  ocean  waters  is  dominated 
by  two  processes:  absorption  and  scattering  by  nonabsorbing 
particles  and  plankton.  Measurements  of  the  differential  scattering 
cross  section  show  that  the  relative  angular  dependence  of  the 
scattering  cross-section  is  remarkably  constant  even  though  the 
total  volumetric  scattering  cross  section  varies 

substantially.  This  suggests  that  we  represent  the  volumetric 
differential  scattering  cross  section  B(z,0)  by  the  separable  form 


e(z,0)  -  0T(«)  X(0)  (1) 

where  z  denotes  depth.  Figure  1  shows  that  the  observed  angular 
dependence  of  the  scattering  can  be  well  represented  by  the  formula 


0n 

X(0) - 2 -  (2) 


0q  ■  .12  radians  (7  degrees) 


where  X(0)  satisfies  the  small-angle  normalization  condition 


Figure  1.  Measurements  of  volumetric  differential  scattering  cross-section  and  the  theoretical 
fit,  formula  (2).  [From  Fig.  IS,  N.  G.  Jerlov,  Marine  Optics,  Elsevier,  (1976),  p.38.] 


(3) 


x(e)2ir0de 


1  . 


Formula  (2)  fits  the  data  well  in  the  region  0  <  0  <  i  ;  the 

2 

contribution  from  angles  in  the  region  —  <  0  <  it  to  the  total 

2 

scattering  cross  section  is  small* 


Let  us  note  that,  while  the  normalization  Integral  (3) 
converges  nicely,  the  mean-square  scattering  angle 


02  > 


(4) 


is  divergent  when  small-angle  approximations  are  used.  For  this 

1-3 

reason,  the  previous  treatments  which  were  based  on  a  small  mean- 
square  scattering  angle  and  a  Fokker-Planck  diffusion  in  angles  are 

not  appropriate.  In  this  work,  we  do  not  use  a  Fokker-Planck  model, 
but  retain  the  small  angle  approximation  to  obtain  a  direct  solution 
to  the  radiative  transfer  equation.  A  closed  expression  for  the 
rate  of  radial  spreading  of  an  initial  pencil  beam  is  found.  For 
laser-radar  applications,  the  scaling  of  the  range  straggling  of  a 
pulse  can  be  found  analytically,  but  the  actual  distribution  comes 
from  a  numerical  solution  to  a  partial  differential  equation. 


RADIATIVE  TRANSFER  EQUATION 


The  radiative  transfer  equation  in  the  small-angle 
approximation  can  be  written 


“li  +  li+e  .  -  -(cH-RpJf 

C  at  3z 

+  sT  JdV  f(e')  x(e  -  e'; 


where  a(z)  is  the  absorption  coefficient,  &T(z)  is  the  total 
scattering  cross  section,  and  n  is  the  index  of  refraction.  The 
variation  of  a  and  3  with  depth  z  is  arbitrary.  The  solution 


to  (5)  takes  the  form 


g|-  t-zj  ♦(z,r,©) 


where  the  time  dependence  of  g  can  be  arbitrary  and  ^  satisfies 
the  time-independent  equation. 


11  +  e  .  -1  i>  -  -(or*-eT)^  +  eT  f dV*(z,r,e')  x(e-e')  (7) 
3z  ,♦  J 


Equation  (7)  can  be  solved  by  Fourier  transform.  Introducing  the 
definition 
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*(z,;,e) .  f j*L  f  Ji.  eip-®  ft.,!,;) 

J  Oif'i2  J  r>i r'*2 


one  arrives  at  the  equation 


- k  •  il  -  -(crfft.)  *  +  eT  ¥  X(p) 


The  Fourier  transform  of  the  scattering  function  depends  only  on  the 
magnitude  of  p  * 


-  f d2e  «-lp'9 


x(e)  -  2ir  J  ©de  •  JQ(pe)  x(e) 


I  f  Jo(po0u)  -  i  / 
T  i)  i+u2 


-Penv 


0  (i-v2) 


Sections  6.532,  8.551,  and  8.43  of  Gradshteyn  and  Ryzhik^  were  used 
to  transform  the  integral. 


Equation  (9)  is  hyperbolic  ahd  is  solved  by  integrating 
along  characteristics 


?(z,k,p)  -  exp  -tq  +  f  8T(z')[x(s)-l]dz' 


where 


T 

a 


a(z' )dz 


•¥  *► 

s  ■  |p  +  k  (z-z' ) | 


(12) 


(13) 


when  z  *  0  ,  ¥  takes  on  the  value  ¥  *  1  corresponding  to  a 

delta-function  beam  i|>  ■  6(r)  S(®)  •  Equations  (8),  (10),  (11)-(13) 
constitute  a  solution  of  the  small-angle  radiative-transfer  problem. 


APPLICATIONS 


In  general,  an  evaluation  of  the  integrals  (8)  and  (11) 
calls  for  numerical  work.  But  fortunately,  two  Important  special 
cases  can  be  evaluated  analytically.  Our  interest  centers  on  the 
case  where  the  optical  depth  to  small-angle  scattering 

=  f  BT(z')dz'  (14) 

*0 

is  large.  In  this  limit,  we  can  utilize  the  Taylor  expansion  of 
X(s) 

X(s)  -  1  -  1  80  +  i.  g202  +  .  .  .  .  (15) 

n  U  4  U 

obtainable  directly  from  the  last  representation  of  equation  (10). 
Only  the  first  term  of  this  series  need  be  retained,  and  we  find 
that 


'f(z,k,p)  -  exp 


BT(z')dz'  |p  +  k  (z-z')| 


(16) 


The  first  special  case  concerns  the  intensity  as  a  function 
of  radius  integrated  over  all  angles 


9 


♦  f  ik*r..  ♦  ♦  ♦ 

I(z»r)  “  /  d0  ♦(Z»r>0)  "  /  * - —  *  (z>k»P) 


(2w) 


p-0 


kdk  jQ<kr) 


2* 


exp 


(‘  7  kV ) 


(17) 


where 


z 


/*  0T(z')(z-z')dz' 
0 


T 

^  drg'  [*-*(tb')]  . 


(18) 


In  the  special  case  where  0^  Is  uniform,  one  finds 


z 


2 

z 


(19) 


The  Integral  in  (17)  can  be  carried  out  [Gradshteyn  and  Ryzhik*, 
section  6.623]  to  give 


l(z*r)  - 


(20) 


Equation  (20)  governs  the  spreading  of  an  initial  pencil  beam  as  a 
result  of  multiple  scattering.  The  radius  increases  as  the  square 
of  the  depth,  rather  than  the  3/2-power  of  depth  dependence 


10 


t 


characteristic  of  theories  based  on  a  small  mean  square  scattering 
angle. 


The  second  case  concerns  the  degradation  of  the  image  of  a 
point-source  object  emitting  almost  Isotropically  in  angle.  At  the 
source,  the  intensity  distribution  is  given  by 


♦  ♦  -e2/2(Ae)2 

*(r,e)  -  6(r)  5 - 

2*(A0)2 


The  angle  A0  has  no  intrinsic  significance;  it  is  chosen  to 
satisfy  0q  <  <  A0  <  1  so  that  the  small-angle  approximation 
remains  valid  while  the  initial  distribution  is  diffuse  compared  to 
the  typical  scattering  angle  0Q  .  Direct  integration  yields  the 
initial  Fourier  transform 


_2p2 

P  eo 


*(k.p)  - 


Starting  with  the  Initial  condition  (22),  one  can  solve  (9)  by 
integrating  along  characteristics  to  obtain 


2  2 

-  ♦  s n  U0> 

*(*>fc,p)  “  exp  -tq - - - 


+  jZ  ST(z')[x(s)-l]dz' 
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sQ-  IP  +  kz| 


(24 


i 


Figure  2. 

f 


Geometry  relevant  to  degradation  of  a  point-source  image.  The  origin  of  z  is  at 
the  point  source.  The  vector  0  is  the  component  of  a  unit  n  vector  which  is 
perpendicular  to  the  z-axis.  The  imaging  system  is  far  removed  and  records  only 
radiation  with  (T  =  0. 


To  Imaging 
System 


l 


At  this  point,  it  is  useful  to  regard  A0  as  formally  large  which 
permits  one  to  ignore  p'  relative  to  kz'  in  the  integral 


involving  $j(z')  .  Thus  one  obtains 


♦  ♦  Ta  f  , 2 ?  20nk  rz 

♦(z,r,e)  - — - -  I  -  exp  ik«(r-ze)  -  i  &_(z')z'dz' 

2x(A0)2  J  (2w)2  *  0 


where 


2’<4e)2  «*[|;  -  ze|2+  (4  ;  e0)2f 


-  C  eT(z')z'dz'  . 


Note  definitions  (18)  and  (29)  for  z  are  identical  when  the 
different  origin  of  z  is  taken  into  account*  The  factor 
l/[2w(A0)2]  represents  simply  the  intensity-per-steradian  of  the 
point  source* 


The  image  formed  according  to  the  system  portrayed  in 

♦ 

Figure  2  depends  on  an  evaluation  of  (28)  at  0  -  0  *  Hence, 
multiple  scattering  degrades  the  image  of  a  point-source  into  one 
with  the  finite  size  Ar  -  (2z0q/it)  . 
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It  is  important  to  emphasize  that  our  approximate  methods 
are  appropriate  to  the  large  optical  depth  regime 

t  -  fZ  e  (z')dz'  >  1  .  (30) 

*/o  1 

But  even  in  this  regime,  there  will  be  an  unscattered  component  of 
strength  e~T  which  could  appear  much  brighter  in  an  imaging  system 
if  there  is  a  true  point  source.  In  experiments  dealing  with  the 
propagation  of  laser  beams,  this  bright  core  dominates  the 
experimental  results^ >^,  even  though  the  fraction  of  photons  in 
the  core  decreases  exponentially.  But  for  extended  objects  of 
finite  surface  brightness,  the  image  will  be  dominated  by  the 
multiplely-scattered  component  which  contains  most  of  the  intensity. 


RANGE-STRAGGLING  OF  A  LASER  RADAR 


The  range  resolution  of  a  laser  radar  in  a  multiple- 
scattering  medium  is  degraded  because  the  photons  do  not  travel  in 
straight  lines.  Let  us  write  the  distribution  function  f  as 


f  -  f 


(31) 


*2 


(32) 


and  include  second-order  multiple  scattering  effects  in  the 
radiative  transfer  equation 


c  3t  \  2/  3z  +• 


3r 


/«  ♦  ♦  ♦  ♦  ♦ 

dz  e'  f(z,*2,r,e')  x(e-e') 


(33) 


Range-straggling  appears  through  the  0^/2  —  term.  In  any  radar 

3s 

system. 


11-  >  >  11  (34) 

3*2  3* 


and  multiple-scale  arguments  permit  one  to  simplify  (33)  to 


i£  +  .£j£+e.M  -  -(cri-eJf  +  eT  fd2  e'f(e')  x(e-e) 
3z  2  3z,  ^  J 


Analytic  progress  can  be  made  if  we  restrict  our  attention  to 


g(z,z_,e)  defined  by 


g(z,*2,e)  •  f d2  r  f  .  (36) 

Physically,  %2  represents  the  distribution  function  governing  a 
plane  wave  propagating  vertically  downward.  Clearly  g  satisfies 
the  equation 

M  +  ®i  i* - +  r.  fi2  e-  g(e')  x(e-e')  (37) 

3z  2  3z2  J 

which  can  be  Fourler-transforned  in  angle  to  be 


l£-i_i_pi_i£_  -  gTc(p) [x(P)-i ] 

3z  2p  3p  dp  3z2 

where 


/♦  «► 

de  e-ip*e  g(z,z2,e) 


The  Taylor  expansion  of  X(p)  results  in 
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(40) 


I 

I 


! 

S 


—  -  —  —  p  —  —  -  -6-  -  e.pG  • 

3z  2p  dp  3p  3z2  ir 

Our  delta-function  Initial  conditions  permit  the  use  of  similarity 
arguments  when  is  constant.  Thus,  we  can  Introduce  the 

similarity  variables 


u  -  pBT20oz/ir 

9  l  _1 

v  -  2z2[(BT20o/ir)/  zJ] 

F  -  G[(0T2OO/*)2  z3]  (41) 

and  transform  (40)  into 

3  ±_  (VF)  +  I  i-  u  1-  -  uF  +  u  *1  •  (42) 

3v  u  3u  3u  \3v/  3u 

The  function  F  will  vanish  for  v  <  0  and  as  v  •*  .  An 
integration  with  respect  to  v  yields 


(43) 

(44) 


Along  the  line  v  ■  0  ,  the  solution  Is  F  ■  constant  .  The 
constant  will  be  determined  by  the  normalization  condition 


f  Fdv  -  e"11 

which  is  readily  derivable  from  (42). 


0 


Further  progress  in  solving  (43  -  44)  does  not  appear 
possible  without  computational  work.  But  the  similarity  scaling 
(41)  shows  that  the  range  resolution  scales  as 


An  accurate  computational  solution  of  (42)  depends  on 

obtaining  the  correct  starting  values  as  v  ♦  0  and  as  v  *  •>  . 
Near  v  »  0  ,  we  can  use  the  trial  solution 


F  "  v»  Ft(y) 


y  ■  u  v^/^  . 


(< 


Substitution  of  this  form  into  equation  (42)  yields  the  result  as 
v  ♦  0 


y  3y  3y 


y  Fi 


where  the  correction  terms  are  of  order  v*'3.  por  iarge  y  ,  the 
solution  takes  the  form 


F 


1 


A  sin 


(49) 


when 


n  - 


In  the  limit  y  ♦  0  ,  the  desired  solution  has  the  form 


F 


1 


1  + 


9(n+l) 


(50) 


These  three  free  parameters  A,  B,  n  permit  one  to  solve  for  a  , 
and  to  liminate  the  two  divergent  solutions  near  the  origin. 


In  the  limit  v  ,  equation  (42)  is  solved  by 


3  — u 
v  e 


(51) 


(vu3) 


where  we  require  s  >  1  in  order  for  the  normalization  condition  (45) 


to  converge.  Let  us  note  that  (51)  is  not  valid  for  u  <  — i —  •  We 

-  ~  yl/3 

also  require  that  /  Fu  du  <  »  ,  leading  to  s  <  5/3  . 


leading  to  s  <  5/3 


SUMMARY 


A  simple  fit  to  the  shape  of  the  differential  cross  section 
allows  substantial  analytic  progress  to  be  made,  including  equation 
(20)  for  the  radial  spreading  of  underwater  intensity,  equation  (28) 
for  the  angle~radius  relationship  for  an  underwater  point  source  of 
illumination,  and  estimate  (46)  for  the  range  resolution  of  a  laser 
radar . 
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